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ABSTRACT 

We investigate the effects of a coupled Dark Energy (cDE) scalar field on the 
alignment between satellites and matter distributions in galaxy clusters. Using 
high-resolution N-body simulations for ACDM and cDE cosmological models, we 
compute the probability density distribution for the alignment angle between 
the satellite galaxies and underlying matter distributions, finding a difference be- 
tween the two scenarios. With respect to ACDM, in cDE cosmologies the satel- 
lite galaxies are less preferentially located along the major axis of the matter 
distribution, possibly reducing the tension with observational data. A physical 
explanation is that the coupling between dark matter and dark energy acts as an 
additional tidal force on the satellite galaxies diminishing the alignments between 
their distribution and the matter one. Through a Wald test based on the gen- 
eralized x 2 statistics, the null hypothesis that the two probability distributions 
come from the same parent population is rejected at the 99% confidence level. It 
is concluded that the galaxy-matter alignment in clusters may provide a unique 
probe of dark sector interactions as well as the nature of dark energy. 

Subject headings: cosmology:theory — large-scale structure of Universe 

1. INTRODUCTION 

The concordance cosmological model - based on a family of non-relativistic massive 
particles possibly weakly interacting with the standard sector of particle physics (Cold Dark 
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Matter particles, CDM hereafter) and on a cosmological constant A - very successfully 
accounts for a large wealth of observational data, but poses significant theoretical puzzles and 
requires an extreme fine-tuning of its basic parameters. In this context, alternative scenarios 
have been explored, as i.e. the possibility that the cosmological constant A be replaced by 
a Dark Energy (DE) compone nt represented by a classical dynamical scalar field <fi evolving 



in a self interaction potential (jWetterichlll988l ; iRatra h Peebles! 1 19881 ). A further step in the 



exploration of alternative cosmological models has been recently proposed, specula ting about 



1995; 


Amendola 


2000; 


Amendola 


2004) 



been widely investigated in recent y ears concerning their background evolution and their 



effects on structure formation (e.g. bylM ainin i fc Bonomettol ; iPettorino fc Baccigalupill2008 
Abdalla et al. 2009 ; Wintergerst fc Pettorino 201 ol ) by means of specifically designed N-body 



(Maccio et al. 


2004; 


Baldi et al. 


2010; 


Baldi 


2010a: 


Li & Barrow 


2010) 



These studies have led to the highlight of some distinctive features in the evolved matter 
density fields of cDE models and have shown how the properties of nonlinear collapsed 
objects can be significantly affected by the d ark interactions in a potentially observable way 



(iBaldi &; Viell I201CH : iBaldi &; Pettorino! 120101 1 . Therefore, one of the most important tasks 
in the present investigation of these alternative models consists in linking the predictions 
obtained for the CDM density distribution to directly observable quantities. 



Led (120101 ) very recently brought up a speculative idea, namely that the dark sector 
interaction might be probed by measuring the alignment between galaxy and matter distri- 
butions in triaxial clusters. In ACDM cosmology, the spatial correlations of the large-scale 
tidal fields yield the strong align ments of the satell ite gal a xies in the clusters with the cluster 
dark matter distributions (e.g.. lAltay et al.l 120061 ) . iLed (120101 ) claimed that the alignment 
of the cluster galaxies with the underlying CDM distribution would be weaker in the cDE 
models than in ACDM since the spatial correlations of the tidal fields inside the clusters 
would be less strong due to the existence of a fifth force. If this is really the case, the align- 
ment between galaxy and matter distributions in triaxial clusters would provide a unique 
probe of the cDE scenarios, especially because the alignment angle is insensitive to the other 
cosmological parameters. 



To back up the speculative idea of ILed (120101 ). it is essential to examine numerically 
whether or not the cDE scalar field truly makes a detectable difference on the cluster galaxy- 
matter alignment. The goal of this Paper is to use high- resolution N-body simulations to 
study how the alignment between satellites and the underlying dark matter distribu tion in 
galaxy clusters changes in cDE models as compared to ACDM (IKomatsu et al.ll2010f ). 
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SIMULATIONS AND DATA ANALYSIS 



The dark interaction between the DE scalar field <fr and CDM which characterizes cDE 
models affects the growth of cosmic structures in a manifold way. First of all, the background 
expansion history of the universe is changed with respect to the standard ACDM evolution 
by the dynamical nature of the DE scalar field <fi. Furthermore, the mass of CDM particles 
changes in time due to the exchange of energy with the DE sector, and the coupling between 
the two fields determines a new long-range fifth force between CDM particles and an effective 
"modified inertia" in the form of a velocity-dependent acceleration. These features can be 
expressed by a modified acceleration equation for CDM particles in the Newtonian limit: 



G[l + 2/3 2 (0)]m i ( ( /»)r 



pi 



where the reduced Planck mass is defined as M P1 = l/y/SnG, /3((p) is the coupling strength, 
Vi is the velocity of the i-th particle, is the vector distance between the i-th and the j-th 
particles, and the sum extends to all the CDM particles in the universe. The first term on 
the right-hand side of Equation [1] represents the "modified inertia", while the [1 + 2 / 3 2 (</>)] 
factor that multiplies the gravitational constant is t he "fifth- force" . All these effects have 
been studied in the literature (see e.g. iBaldil l2010al lbl) - to which we refer the interested 
reader for a more detailed explanation of the cDE effects - and need to be properly included 
in N-body algorithms in order to consistently evolve the growth of density perturbations in 
these scenarios. 

In the present work, we rely on the numeric al implementat ion of cDE m odels in the 



reeP M cosmological N-body code GADGET-2 (jSpringell 120051 ) presented in iBaldi et al. 



(120 lOl ). with which we run a series of cosmological simulations for a fiducial ACDM cos- 



mology with WMAP7 parameters and two cDE scenarios with the same background pa- 
rameters and with linear density perturbations normalized at the present time z — 0. Our 
c DE models consist o f one constant coupling scenario with = 0.2 (the EXP005 model 



of Baldi Sz Pettorinol (hoiO)) and one variable coupling model with /3 = 0.4 ■ e 2< ^/ Mpi (the 



EXP010e2 model of IBaldil (j2010al )). Notice that the coupling definition assumed in this work 
is different from some of the literatu r e and our values of the coupling are a/2/3 smaller than 
the ones quoted in e.g. lAmendola I (120001 ). For these three different cosmologies we run a 
CDM simulation on a cosmological box of 100 comoving /i -1 Mpc aside, with 512 3 particles. 
This corresponds to a space resolution of e s = 4 h^ 1 kpc and to a mass resolution at z = 
of m = 5.6 x 10 8 h~~ l M & . We then link CDM particles into groups with a friends-of-friends 
procedure and for ea ch main halo we iden tify all the bound substructures by means of the 
SUBFIND algorithm ( iSpringel et al.ll200fl ). Let us remind that we use the same initial ran- 
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dom field in all our simulations. This allows us the possibility of a one-to-one comparison 
between different models. 

In order to compare the results of our numerical simulations with observations we need 
to assign to each (sub)halo a stellar mass. For this purpose w e use the recently derived 
stellar-to-halo mass relationship presented in iMoster et al.l (120101 ) . In this paper the authors 
give a series of simple fitting formulae to derive the stellar mass of the galaxy hosted by a 
given dark matter halo as a function of the virial mass (for central galaxies) or as a function 
of the maximum mass over the halo's history (for satellite galaxies). For all our substructures 
we therefore build a full merger tree and we assign to eac h subh alo a stellar mass based on 
the latter criterion by using the formulae of IMoster et al.l ( 120101 ) . 



However, these formulae have been derived by comparing the observed galaxy luminosity 
function with the halo mass function in a ACDM universe. This implies that we cannot 
directly apply them in the case of cDE models. In order to overcome this issue, we first 
find a functio n f(M, z) able to m ap the cDE mass function onto the ACDM one and then 
we apply the IMoster et al.l ( 120101 ) fitting equations. The CDM density and the location 
of the 30 most luminous substructures inside two virial radii for the same massive cluster 
within the three different cosmologies under investigation are plotted in the three panels of 
Figure (U where a stronger tendency of alignment between the CDM distribution and the 
selected substructures in the ACDM model as compared to cDE scenarios is already visible. 



3. ALIGNMENT BETWEEN GALAXY AND MATTER DISTRIBUTIONS 

IN CLUSTERS 

From each cluster catalog obtained from the simulations for the ACDM and cDE cos- 
mologies, we select those cluster halos which have five or more galaxies within some cut-off 
radius r c (< R V i r ) from the halo center. Then, we measure the alignment angles between 
the galaxy and dark matter distributions of each selected cluster halo through the following 
procedure. First, the inertia tensor of the galaxy distribution of each selected cluster halo, 
is evaluated as 

p 

where m* t is the stellar mass of the /3-th galaxy and (xpj) is the position vector of the (3-th 
satellite galaxy relative to the center of the stellar mass of the galaxy distribution. Weighting 
the satellites by stellar mass is important for a comparison with observations since in practice 
the shapes of galaxy distributions are often obtained from the luminosity-weighted satellites. 
Finding the unit eigenvector of (1^) corresponding to the largest eigenvalue, we determine 
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the major axis, e©, of the galaxy distribution of each selected cluster halo. 

Similarly, we evaluate the inertia tensor, (Iff), of the CDM distribution of each selected 
cluster halo as 

lij = ^ ^ Xg^XgJ, (3) 

a 

where (x a ^) is the position vector of the a-th dark matter particle relative to the halo 
center. The summation in Equation Q goes over all dark matter particles located within 
the halo virial radius R V i r . Finding the unit eigenvector of (Iff) corresponding to the largest 
eigenvalue, we determine the major axis, e^-, of the CDM distribution of each selected cluster 
halo. The alignment angle, 9, between dark matter and galaxy distributions of each selected 
halo is now calculated as cos 8 = |e_^f • ec\ where cos 8 is in the range of [0, 1]. Binning 
the values of cos 8 and counting the number of the halos belonging to each cos^-bin, we 
derive the probability density distribution, p(cos8), for the ACDM and cDE cosmologies, 
separately. When the cut-off radius r c for the satellite galaxies is set at 0.8_R„ r (see Figure 
H]), the total number of the selected cluster halos which have five or more galaxies within r c 
is found to be 186 and 116 for the ACDM and cDE cosmologies, respectively. 

Before making a comparison between the probability distributions of cos 8 from the two 
cosmologies, it is worth examining how robust our measurement of cos 8 is. If the dark matter 
distribution of a given halo is spherically symmetric, then the determination of the major 
principal axis of its inertia tensor would be quite ambiguous. Therefore, the robustness of 
the measurement of cos 8 of a given halo critically depends on how spherical its dark matter 
distribution is. The sphericity of a given halo is conventionally defined as S = yu^Jui-, 
where Ui and u% are the largest and the smallest eigenvalues of the inertia tensor of the dark 
matter distribution, respectively. Using this definition, we measure the sphericity of each 
selected halo for both cosmologies. Figure |2] plots the number distribution of the selected 
halos versus their sphericity for cDE and ACDM cosmologies as solid and dashed histogram, 
respectively. As it can be seen, for the case of cDE there is no halo whose sphericity is larger 
than 0.9, while for ACDM there is only one halo whose sphericity exceeds 0.9. Thus, our 
measurement of cos# is robust, free from the ambiguity caused by the spherically symmetric 
dark matter distribution. 

Given that our samples of the halos for both cosmologies are not large enough to be 
fully representative of the true parent populations, we perform the bootstrap error an alysis 



to estimate the uncertainties in the measurement of p(cos8) (IWall fc Jenkins! 120031 ) . We 
first construct 1000 bootstrap resamples and measure p(cos8) from each resample. The 
bootstrap error, <7b, at each bin is then calculated as one standard deviation among the 1000 
remeasurements, = ([p(cosdi) — (^(cos^j))] 2 ) 1 / 2 , where the ensemble average is taken over 
the 1000 resamples and (p(cos#j)) represents the bootstrap mean value at the i-th bin. The 
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top panel of Figure [3] plots the resulting ^(cos^) with the bootstrap errors for cDE and 
ACDM cosmologies as square and triangle dots, respectively. For this plot, the cut-off radius 
r c for the satellite galaxies is set at 0.8-R„j r . 

As shown in Figure [3j the probability density distribution increases with cos 9 in both 
cosmologies, which indicates that the major axes of the galaxy and dark matter distribu- 
tions in clusters tend to be strongly aligned. However, we note a difference between the two 
cosmologies: p(cos8) increases less rapidly in the cDE model. At the fourth cos^-bin (cor- 
responding to 0.6 < cos 9 < 0.8), the difference between the two models exceeds 3cr&. The 
bottom panel of Figure [3] plots the ratio between the two probability density distributions: 
PcDe(cos 6) I Pacdm(cos 6) , which demonstrates clearly that the ratio deviates from unity and 
the degree of the deviation reaches the maximum at the fourth bin. 

Now that we find a difference between the two cosmologies in the tendency of the cluster 
galaxy-matter alignment, we would like to investigate how the result depends on the cut- 
off radius r c . In fact, it is naturally expected that as r c decreases the alignment tendency 
would decrease. The major axes of the dark matter distribution are determined at R V i r 
while the major axes of the galaxy distribution are measured at r c . Therefore, the difference 
Rvir — r c represents the distance by which the inner and outer tidal fields are separated. 
As r c decreases, the spatial correlations between the two tidal fields would decrease, which 
would be manifested by the decrease of the strength of the alignment between the major 
axes of the galaxy and matter distributions in clusters. 

Repeating the whole process described above, we calculate the probability density at 
the fourth bin, p(0.6 < cos8 < 0.8), for five different choices for r c . The top panel of Figure 
S] plots the result for cDE and ACDM cosmologies as square and triangle dots, respectively. 
The bottom panel of Figure @] plots the ratio of the square points (cDE) to the triangle 
points (ACDM). Figure S] reveals that in the cDE scenario the probability density, p(0.6 < 
cos 6 < 0.8), decreases more rapidly as r c decreases. This result may be interpreted as the 
effect of coupled dark energy on the spatial correlations of the tidal fields inside the cluster 
halos. Due to the coupling between dark matter and dark energy, there is an additional 
tidal force exerted on the galaxies inside the clusters, which plays a role of reducing further 
the correlations of galaxy distributions with matter distributions. Furthermore, the velocity- 
dependent acceleration arising as a consequence of the "modified inertia" described by the 
first term on the right-hand side of Equation (JTJ also contributes to an effective modification 
of the internal tidal field of the cluster experienced by the substructures since it determines 
an effective extra force along the direction of motion of the satellite halos. 

To make a direct comparison with observations, it may be more useful to find the prob- 
ability density distribution of the two dimensional projected alignment angles, #2d- For each 
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cosmological model, we project and of each selected halo cluster onto the xy, yz and zx 
planes and calculate the two dimensional alignment angle 9 2 d between the projected major 
axes of dark matter and galaxy distributions for all three projection cases and determine 
the probability density distribution, p{9 2 d), by taking the average over the three cases. To 
estimate the uncertainties in the measurement of p(62d), we also perform the boostrap error 
analysis. The top panel of Figure |5] plots the p{9 2( i) with the bootstrap errors for cDE and 
ACDM cosmologies as square and triangle dots, respectively. As one can see, the two dimen- 
sional projection does not dilute the difference in the alignment tendency between the two 
models. The probability density distribution, p(9 2 d), decreases less rapidly as 9 2 & increases 
in the cDE model than in the ACDM model. Note that for a three dimensional distribution, 
it is p(cos 9) that is expected to be uniformly distributed in case of no alignment while for 
the two dimensional case, it is p(9 2 d). 

To quantify the statistical significance of the difference in the galaxy-matter align- 
ments between ACDM and cDE cosmologies, we perform a Wald test for the ratio, £ = 
PcDe{cos9)/p c deacdm(cos9). In the null hypothesis that the two probability density 
distributions shown in Figures [3] and are in fact from the same parent population, the 
expectation value of £ is unity at all bins. Taking into account the possible cross-correlations 
of p(cos 9) between the different bins, we calculate the generalized x 2 defined in terms of the 
full covariance matrix as 

-I)- ( 4 ) 
Here CZ 1 is the inverse of the covariance matrix, = ((& — — £,)), where the 

ensemble average is taken over the 1000 bootstrap resamples and denotes the mean value 
of £ averaged over the bootstrap resamples at the i-th bin. In Equation (j4j), — 1) expresses 
how much the numerical result of £ at the z-th bin deviates from the expectation value of 
unity. We find \ 2 = 15.08 and 14.97 for the three and two dimensional case, respectively. 
Since the number of degrees of freedom for x 2 is 5, its value corresponds to the null hypothesis 
being rejected at the 98.9% confidence level for both cases. 



DISCUSSION AND CONCLUSIONS 



Using data from high-resolution N-body simulations, we have shown that the ACDM 
and cDE cosmologies yield different strengths of the align ment b etwee n satellite galaxy and 
matter distributions in triaxial clusters. As speculated by iLed (120101 ). it is found that the 
alignment is less strong in cDE cosmologies than in the standard ACDM scenario. The null 
hypothesis that the two cosmologies give the same alignment tendency is rejected at the 
98.9% confidence level through a Wald test based on the generalized \ 2 statistics. 
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The differences in the satellite distributions between ACDM and cDE models may be 
understood as follows. There are two possible me chanisms for expl aining; the alignments 
between the dark matter and galaxy distributions (lAltay et al.l 120061 ) . The first one is the 
spatial correlations between the large scale external tidal fields and the internal tidal fields 
due to the cluster potentials which affect he non-spherical shapes of dark matter compo- 
nents of galaxy clusters and the anisotropic distributions of their satellite galaxies, respec- 
tively. The second one is the filamentary merging/accretion of matter and galaxies along the 
large-scale filaments, which also induce the alignments between the dark matter and galaxy 
distributions in clusters. 

In the ACDM cosmology, the tidal fields evolve nonlinearly only via gravity. Whereas in 
the cDE cosmology the coupling between dark energy and dark matter generate additional 
(more isotropic) tidal effects in the clusters, which weaken the spatial correlations between 
the external and internal tidal fields and also redistribute previously accreted satellites that 
were initially in a flattened distribution. Furthermore, the "modified inertia" characterizing 
cDE models would contribute to an effective enhancement of the internal cluster tidal field. 



Recently lOguri et al.l ( 120101 ) compared the observed distribution of satellite galaxies 
in clusters with the shape of the cluster dark ma tter distribution determined from weak 
lensing, finding a very weak alignment (ILed l2010l ). Due to the small sample (composed 
of only 13 clusters), it is difficult to say at the moment whether or not this observational 
signal supports the cDE scenario. In any case, a crucial implication of our results is that the 
alignment between galaxy and matter distributions in clusters is in principle a unique probe 
of dark sector interactions. 
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Fig. 1. — CDM projected density (grey scale) and location of the 30 most luminous sub- 
structures (colored squares) within i? 2 oo for the most massive cluster formed in our simulated 
ACDM (left), EXP005 (middle) and EXP010e2 (right) cosmologies, respectively 
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Fig. 2. — Number distributions of the selected halos as a function of the sphericity of their 
dark matter distribution for two different cosmologiesxDE (solid histogram) and ACDM 
(dashed histogram). 
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Fig. 3. — Probability density distribution of the cosines of the angles between the major 
axes of dark matter and galaxy distributions of the simulated clusters for two different 
cosmologies: cDE (solid line+square dots) and ACDM (dashed line+triangle dots) in the 
top panel. The errors represent the la scatter among the 1000 bootstrap resamples. The 
ratios of the square points to the triangle points are plotted as solid line in the bottom panel. 
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Fig. 4. — Probability density that the cosines of the angles between the major axes of dark 
matter and galaxy distributions of the simulated clusters are in the range of (0.6, 0.8) with 
the bootstrap errors as a function of the cut-off radius of the satellite galaxies r c for two 
cosmological models: cDE (solid line+ square dots) and ACDM (dashed line+triangle dots), 
in the top panel The ratios of the square points to the triangle points at each bin are plotted 
as solid line in the bottom panel. 
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Fig. 5. — Probability density distribution of the angles between the two dimensional pro- 
jected major axes of dark matter and galaxy distributions with the bootstrap errors for two 
different cosmologies: cDE (solid line+square dots) and ACDM (dashed line+triangle dots). 
The ratios of the square points to the triangle points are plotted as solid line in the bottom 
panel. 



